Bilayer Hubbard model for ^He: a cluster dynamical mean-field calculation 
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Inspired by recent experiments on bilayer '^He, we consider a bilayer Hubbard model on a trian- 
gular lattice. For appropriate model parameters, we observe a band-selective Mott transition at a 
critical chemical potential, pc, corresponding to the solidification of the fermions in the first layer. 
The growth of the effective mass on the metallic side (}i < pc) is cut off by a first order transition 
in which the first layer fermions drop out of the Luttinger volume and their spin degrees of freedom 
become locked in a spin singlet state. These results are obtained from a cluster dynamical mean-field 
calculation on an eight-site cluster with a quantum Monte Carlo cluster solver. 



The solidification of '^He monolayers [1] has been in- 
terpreted as a density-driven Mott transition in which 
the effective mass diverges [21 E] . Below the critical den- 
sity, the metallic phase corresponds to a nearly localized 
Fermi liquid and, beyond the critical density, to a spin- 
disordered solid. Recently, it has been possible to realize 
bilayers of '^He [1] (atop a frozen ^He substrate, itself 
adsorbed onto graphite) with the special property that 
the second layer begins to form before the first has solidi- 
fied. Since the first layer is close to a Mott transition, the 
•^He fermions in this layer are slow (i.e. heavy), whereas 
those in the second layer are fast. This combination of 
fast and slow dynamics — corresponding to wide and nar- 
row fermion conduction bands — is completely analogous 
to the situation in electronic heavy fermion materials, al- 
beit without the complication of crystal field and spin 
orbit effects. 

According to this picture, prior to solidification of the 
first layer, one expects an enhanced effective mass and a 
Luttinger volume that counts both the first- and second- 
layer populations. Moreover, one naively anticipates that 
further '^He deposition will eventually cause the effective 
mass to diverge. In experiment, the effective mass is 
indeed observed to increase as a function of the total 
•^He concentration, but its growth is interrupted by an 
intervening phase H] in which the first "^He layer is a spin- 
disordered insulator, decoupled from the second layer. 

The motivation for this Letter is to consider a simple 
lattice model that goes a good way toward reproducing 
the aforementioned experimental situation. As shown in 
Figs, [ija) andjljb), we adopt a stacking of billiard balls 
modeling of bilayer ^He on a triangular lattice defined 
by ai = (l/2,v^/2,0) and = (1,0,0). Each unit cell 
accounts for two ^He positions, Vf = and = |ai — 
|a2 -I- (0,0,03), measured relative to the lattice. Note 
that this presupposes a particular stacking arrangement 
for the second '^He layer. 

Our model can be viewed as a honeycomb lattice whose 
incquivalcnt sites (corresponding to ^He positions in the 
upper and lower layers) are populated by two species 
of fermion, which we label c and /. The tight-binding 
parameters include a nearest-neighbour (interlayer) hop- 
ping V and next-nearest-neighbour (intralayer) hoppings 




FIG. 1: (a),(b) Stacking of billiard balls modeling of bi- 
layer '^He, top and side view, with the *He substrate shown in 
white, (c) Tight-binding modeling with hoppings tc, tf, and 
V. The grey bars connecting / fermion sites (first-layer ^He 
positions) indicate a possible singlet pattern compatible with 
the eight-site 2ai x 2a2 and four-site 2ai x a.2 unit cells, (d) 
The hexagonal Brillouin zone of the triangular lattice and the 
rectangular zone that results from the folding bi |bi- 

tc and tf. See Fig. [ijc). With the inclusion of onsite 
Coulomb repulsion terms, the Hamiltonian reads 
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Here, the mixing element ^(k) = V^(l-l-7k)"'^^^ and the 
dispersion £c(k) = — ic7k + £c expressed in terms 
of the connection 7k = cos(k • ai) -I- cos(k • a2) -I- cos[k • 
(a2— ai)] of the underlying Bravais lattice. The operator 
fic.i ~ X^o- '^i o-'^i o' ^^"^ local '^He density in the upper 
layer. Similar definitions hold for e/(k) and nf,i. 

Except for the complication of the layer stacking (and 
the resulting k-dependent hybridization), this bilayer 
Hubbard model reduces to the Periodic Anderson Model 
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FIG. 2: Single particle occupation number for the / and c 
fermions as a function of temperature. 

as tf 0, a limit in which the bare mass of the / 
fermions diverges. Similar models (though ones with 
a somewhat unrealistic direct layer stacking) have been 
considered for the description of bilayer ■^He in Refs. 
and E] within a slave boson mean-field calculation. Here 
we go one step further and perform calculations within 
the cellular dynamical mean field theory (CDMFT) [7] 
approximation using a supercell defined by the lattice 
constants LcSli and LcB.2 with Lc = 2. Since the original 
unit cell contains two orbitals, this amounts to a total 
of eight orbitals per supercell. The resulting single par- 
ticle Green function, G{'K,iuJm), is a 2L^ x 2L^ matrix 
with crystal momentum K in the Brillouin zone of the 
supercell lattice. The CDMFT calculation involves ne- 
glecting momentum conservation and thereby obtaining 
a K-indepcndcnt self-energy Il(iujm)- This quantity is 
extracted from an Lc x Lc cluster of unit cells embed- 
ded in a dynamical mean field that is determined self- 
consistently. We have solved this cluster problem using 
a standard Hirsch-Fye approach and have symmetrized 
the cluster Green function to obtain the corresponding 
quantity on the lattice: 

" a.l3 

(2) 

Here Xq, denotes the unit cell positions within the super- 
cell, and 1/ run over the c and / orbitals within each 
unit cell, and k and K differ by a reciprocal lattice vec- 
tor of the supercell Bravais lattice. The rotation to real 
frequencies was carried out with a stochastic analytical 
continuation technique [8, 9J. 

We consider the following model parameters: tc ^ tf = 
t, Uc/t =Uf/t= 12, Vo/t = 1/2, e"c/t = 3 and ej/t = 0. 
We have chosen large values of Uc and Uf to reflect the 
contact repulsion of the ^He atoms and to guarantee that 
each single layer is well within the Mott insulating phase 
at half-band filling [ID]. These values of the Hubbard 
interaction lead to low double occupancy, thus gener- 
ating local moments. The difference e° — e° > is a 
crude accounting for the van der Waals forces (both ^He- 



•^He and '^He-'^He) that preferentially fill the first layer. 
Fig. |2] plots the layer densities {hj) and (he) as a func- 
tion of the chemical potential, which controls the overall 
■^He concentration. Analysis of the temperature depen- 
dence of (hf) is consistent with a zero-temperature jump 
of this quantity ai — ^c — — 1.8i. In contrast, {fic) 
increases smoothly with the chemical potential. Since 
^ = {uc + Uf), the jump in the total fermionic density 
signals a density-driven first order transition. In a canon- 
ical ensemble, states with total density lying within the 
jump are phase separated. The first order nature of this 
phase transition can be confirmed explicitly on smaller 
four-site clusters, which can be simulated at much lower 
temperatures before the negative sign problem becomes 
unmanageable [llj . 

The nature of the distinct metallic phases on either 
side of the transition is best understood in terms of the 
single particle spectral function, 

A(k,cj) = -ImTrG(k,w), (3) 

plotted in Fig. [3j For fi < fic, and as exemplified by 
the data set at ^ = —2.5t, the low-energy coherent fea- 
tures of the spectral function compare favorably with a 
slave boson approximation leading to mass-renormalized, 
hybridized bands. This state has a Luttinger volume 
that includes both / and c fermions, and the band with 
the largest Fermi volume has dominant / character. As 
a function of the chemical potential, the effective mass 
of the / band grows; beyond fic, the / band drops 
out of the low-energy physics altogether. This can be 
understood at the static mean field level by a conven- 
tional slave boson theory in competition with local sin- 
glet formation in the first layer. The transition is sig- 
naled by the appearance of an anomalous expectation 

value Ay ~ itj/Uf)J2aiflafj,'^)- The upper inset of 
Fig. [3] shows the band structure that results when this 
singlet order parameter has the configuration depicted in 
Fig. [Tic). 

We understand this transition to be of the band- 
selective Mott type, in which a half-filled band with dom- 
inant / and sub-dominant c character is gapped beyond 
fic- This interpretation is supported by the fact that, be- 
yond He, (fif) does not saturate to unity, as in the case 
of an orbital-selective Mott transition |T5]. The data in 
Fig.[3]show the typical Mott-Hubbard transfer of spectral 
weight from the upper band a,t to ^ 6t = —l.75t) down 
to the Fermi energy as a function of decreasing chemical 
potential [T^]. In the generic Mott-Hubbard scenario, 
doping occurs when the chemical potential reaches the 
lower Hubbard band, which in the present case is situ- 
ated at roughly oj ~ — 6< for fi = — 1.75t. In the band- 
selective Mott transition, doping is provided by changing 
the occupation of the ungapped band. 

To accurately estimate the effective mass renormaliza- 
tion at /i < He, we consider the self-energy, S//(k, ujm), as 
defined by Gj^Ck^icOm) = iw™ - £/(k) H- /i - E//(k, cj,„), 
with the Green function taken from Eq. (^. Since at 
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fjL < He the self-energy is dominated by its frequency de- 
pendence, the effective mass renormalization is very well 
approximated by the inverse quasiparticle residue, 



z-i(k) = i- 



(4) 



This quantity is plotted in Fig. [4] for Fermi wave vec- 
tors of the / band along the F-M and K-F directions in 
the Brillouin zone. A monotonic increase of the effective 
mass is observed, but its divergence is preempted by the 
first order transition ai fi — fic- 

The growth of the effective mass corresponds to a de- 
crease of the coherence temperature, Tcoh, below which 
Fermi liquid behavior manifests itself. To illustrate this. 



A{k,uj), /3t ==20, /i= -1.75t 




A(k,uj), pt = 20, n= -1.875t 
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FIG. 3: The amplitude of the single particle spectral func- 
tion, as defined by Eq. (|3|, is plotted for various values of 
the chemical potential. The lower inset shows the mean-field 
band structure consisting of two quasiparticle bands of mixed 
c and / character. The upper inset shows the single c-only 
band completely decoupled from the gapped, nearly fiat band 
of the singlet-bound / fermions. 




FIG. 4: Quasiparticle residue as defined by Eq. (j4]| and the 
spin scale Ts as a function of chemical potential. The data is 
extracted from simulations at A^c = 8 and I3t — 20. 



we have computed the local spin susceptibility on the 
cluster, as defined by 
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dre*^^'"-(SHx,r).SKx)). (5) 



Below the coherence temperature, x/(*^m = 0) is ex- 
pected to be temperature independent. On the other 
hand, in the temperature region Tcoh < T <C it 
should exhibit Curie- Weiss behavior — the signature of lo- 
cal moment formation. Precisely this behavior is seen in 
Fig. |5] As the chemical potential grows from /i = —2.5t 
to H = — 1.875i, the crossover temperature scale be- 
tween the Curie- Weiss-like and temperature-independent 
Xf{i^7n — 0) tracks the decrease of the inverse effec- 
tive mass and coherence temperature. Beyond the phase 
transition, fi = —1.75t > fj.c, only Curie-Weiss behavior 
is apparent in the considered temperature range. 

Fig. |6] shows that the band-selective Mott transition is 
linked to a sudden growth of the antiferromagnetic cor- 
relations between nearest-neighbor / fermions and to a 
decrease in the intracell c-f spin-spin correlations. This 
data supports the picture that, in the band-selective 
Mott insulating phase, the / quasiparticles are bound 
into spin singlets amongst themselves. The gapping of 
the spin and charge degrees of freedom of the / quasi- 
particles at /i > /ic allows for the decoupling of / and c 
quasiparticles: a c quasiparticle at the Fermi level cannot 
scatter off an / quasiparticle due to the absence of phase 
space. 

The local dynamical spin structure factor. 



bf(u;) = Im 



(6) 



(see Fig. [6]) shows a depletion of spectral weight at low 
energies on both sides of the transition and a consider- 
able sharpening of the line shape in the band-selective 
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FIG. 5: Inverse local spin susceptibility as a function of tem- 
perature for the Nc = 8 cluster. 




FIG. 6: (a) Intracell spin-spin correlations between c and / 
fermions. (b) Nearest-neighbor spin-spin correlations between 
/ fermions. (c) Local dynamical spin structure factor. 



Mott insulating state. At /i < /ic, we can interpret the 



data within an itinerant fermion picture where the mass 
enhancement prior to the band-selective Mott transition 
is taken into account by a renormalization of the hy- 
bridization V and hopping as in a slave boson approach 
[S]. Following this modeling, the peak position in Sf{uj), 
which we will denote by Ts, is expected to track the 
coherence temperature or, equivalently, the inverse effec- 
tive mass. The quantity Ts is plotted in Fig. |4] alongside 
Z(k/) and confirms the above expectations to a good de- 
gree. At n > /Jc, Sf{uj) should be interpreted within a 
localized / fermion picture, in which case the peak po- 
sition is a measure of the excitation energy required to 
break the singlet state of the / fermions. 

We can summarize our results using the terminology 
of heavy fermions [T3]. The band-selective Mott transi- 
tion corresponds to a Kondo breakdown in which the / 
fermions drop out of the Luttinger volume. Within our 
model, it appears that this transition is first order: the 
reduction of the coherence temperature and the enhance- 
ment of the effective mass is interrupted by the formation 
of a spin gapped Mott insulating state of the / fermions. 
The nature of this band-selective Mott insulating state 
is very dependent on the cluster topology. In this work, 
we have considered only clusters with an even number 
of unit cells — thereby implicitly allowing for spin gapped 
insulating states of the / fermions within the CDMFT 
calculation. Despite the breaking of translation invari- 
ance inherent to the CDMFT, the Luttinger sum rule 
still holds when formulated in the Brillouin zone of the 
supercell Bravais lattice. We can only speculate as to the 
nature of this state when the cluster size diverges, but we 
cannot exclude the intriguing possibility that it smoothly 
connects to fractionalized Fermi liquids [T3] [T^ with no 
lattice and spin symmetry breaking but a Luttinger vol- 
ume encompassing only the c fermions. 
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